OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



First Evidence for the Presence of Iron Oxidizing 
Zetaproteobacteria at the Levantine Continental Margins 

Maxim Rubin-Blum 1 *, Gilad Antler 2 , Rami Tsadok 1 , Eli Shemesh 1 , James A. Austin, Jr. 3 , 
Dwight F. Coleman 4 , Beverly N. Goodman-Tchernov 1 , Zvi Ben-Avraham 1 ' 5 , Dan Tchernov 1 

1 The Leon H. Charney School of Marine Sciences, University of Haifa, Haifa, Israel, 2 Department of Earth Sciences, University of Cambridge, Cambridge, United Kingdom, 
3 Institute for Geophysics, Jackson School of Geosciences, University of Texas at Austin, Austin, Texas, United States of America, 4 Graduate School of Oceanography, The 
University of Rhode Island, Narragansett, Rhode Island, United States of America, 5 Department of Geophysical, Atmospheric and Planetary Sciences, Faculty of Exact 
Sciences, Tel Aviv University, Ramat Aviv, Israel 



Abstract 

During the 2010-2011 £/V Nautilus exploration of the Levantine basin's sediments at the depth of 300-1300 m, densely 
patched orange-yellow flocculent mats were observed at various locations along the continental margin of Israel. Cores 
from the mat and the control locations were collected by remotely operated vehicle system (ROV) operated by the £/V 
Nautilus team. Microscopic observation and phylogenetic analysis of microbial 16S and 23S rRNA gene sequences indicated 
the presence of zetaproteobacterial stalk forming Mariprofundus spp. - like prokaryotes in the mats. Bacterial tag-encoded 
FLX amplicon pyrosequencing determined that zetaproteobacterial populations were a dominant fraction of microbial 
community in the biofilm. We show for the first time that zetaproteobacterial may thrive at the continental margins, 
regardless of crustal iron supply, indicating significant fluxes of ferrous iron to the sediment-water interface. In light of this 
discovery, we discuss the potential bioavailability of sediment-water interface iron for organisms in the overlying water 
column. 
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Introduction 

The recently described class of marine microbial iron oxidizers 
(FeOB), the zetaproteobacteria [1], is represented by a stalk- 
forming prokaryote Mariprofundus ferrooxydans [2] , prevails in the 
marine iron-oxidizing biofilms [1,3-7], and is mainly linked to 
hydrothermal activity [4,5,8] and exposed oceanic crust [9,10]. 
Moreover, zetaproteobacteria were recendy shown to participate 
in corrosion processes in the nearshore environments [3,6] as well 
as exist in saline estuaries [1 1] . The autotrophic zetaproteobac- 
teria may play a considerable role in global Fe and carbon cycling 
[4], although their biogeography and metabolic potential are still 
not fully understood. Since zetaproteobacteria studies are mostly 
limited to the deep sea, the iron-dependent carbon fixation is yet 
to be suggested as a considerable source of organic material in the 
benthic environment, and therefore the effects of FeOB on 
benthic-pelagic iron exchange are yet unexplored. 

At sites of hydrothermal fluid emission and exposed ocean crust, 
high fluxes of ferrous iron from the crustal source provide 
metabolic energy to FeOB [4]. In the absence of the latter, 
diagenetic chemical reactions at the sediment-water interface and 
at the top sediment layers control fluxes of important metabolites, 
such as reduced iron [12,13]. Recently, significant benthic iron 
fluxes were recognized in the iron isotopic values measured within 
in situ benthic chambers, which emphasized the importance of the 
microbial iron reduction within sediment for the iron supply to the 



sediment surface and the water column [12]. Given significant 
benthic iron fluxes and oxygenated conditions at the sediment- 
water interface, a niche for the FeOB can be formed at the 
continental margin, regardless of the crustal iron source. 

The net outward flux of soluble iron from the sediment to the 
upper water column layers under upwelling or mixing conditions, 
measured at the continental margins, can provide supplementary 
nutrition to the primary producers and diazotrophs [12,14—16]. 
These findings challenge the paradigm stating that the iron supply 
to the surface waters is controlled by the iron-rich dust supply 
[12,17,18]. In turn, the flux of soluble iron is controlled by the 
organic load that is causing rapid depletion of oxygen and switch 
to manganous and ferruginous respiration closer to the sediment- 
water interface, or by bioturbation/bioirrigation that allow for a 
rapid advection of the ferrous iron to the sediment-water interface 
[12,13]. Moreover, physical or biological resuspension of the 
sediment also leads to the release of soluble particles to the 
medium overlaying the sediment [13,19]. On the other hand, the 
microbially-enhanced iron oxide precipitation at the sediment- 
water interface can alter the bio-availability of the sediment- 
recycled iron, due to different solubility properties of various iron 
oxides [20]. FeOB precipitate nanoparticulate ferrihydrite - like 
phases with short-range structural order [21], although more 
organized iron minerals, such as lepidocrocite, were found in 
association with FeOB [22]. FeOB formed iron oxides are known 
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to be excellent substrates for the iron-reducing bacteria [9], yet 
their bioavailability to other species is still unknown. 

During the 2010-2011 exploration season of the E/V Nautilus 
we studied the deep benthic environment of the Levantine basin, 
the most oligotrophic part of the Mediterranean Sea [23]. The 
contrast between the high iron demand for required for nitrogen 
fixation [24] and the low iron availability in the surface waters [25] 
may explain the low nitrogen fixation rates in this area [26]. The 
paradigm of Fe-rich dust supply to surface waters as the dominant 
iron source includes Levantine basin [27]. In this study we show 
that zetaproteobacterial mats are present throughout Israel's 
continental margins, at depths of 300-1000 m and bring evidence 
showing that such mats are widespread along the Mediterranean 
continental margins. Our findings hint at the potential of the 
FeOB to alter sediment iron bioavailability to the water column. 
Moreover, we provide indirect evidence for the presence of 
significant ferrous iron fluxes from the sediment to the sediment- 
water interface at the marginal areas. 

Materials and Methods 

No specific permissions were required for the locations/ 
activities used in this study. We confirm that the field studies did 
not involve endangered or protected species. The samples were 
collected from an unprotected area: GPS coordinates defining the 
study area are: 32°44. 0945' North, 34°47. 0162' East; 32°42. 
2631' North, 34°33. 5556' East; 32°35. 3227' North, 34°42. 
9271' East; 32°34. 7468' North, 34°44. 4832' East. 

ROV imaging and sample collection 

Nautilus E/V is equipped with Hercules and Argus Remotely 
Operated Vehicle (ROV) systems, which are able to collect high- 
resolution video, oceanographic data, and precision sampling. The 
yellow biofilms at the Levantine sea floor were recorded during 
two legs of the Nautilus E/V, in September 2010 and November 
2011 (NA-009 and NA-019) using the high definition imaging 
system mounted on the Hercules ROV. The areas explored were 
the Achziv canyon at the depths of 500-1 100 m; two locations in 
the Acre area, one deep (1200-1700 m) and one shallower (1000- 
1 200 m); Dor disturbance (both legs, 300-900 m) and Palmachim 
disturbance (both legs, 600-1300 m) [28,29] (Fig. 1). Each 
contiguous survey covered approximately 5 km distance (e.g. 
Fig. lb), while the field of view diameter was ca. 4 m, yielding ca. 
20000 m 2 transect area. Bulk sediment (approximately top 5 cm) 
samples were collected during 20 10 Nautilus E/V field season from 
the Dor disturbance at 32°44. 0945' North, 34°47. 0162' East at 
557 m and from a biofilm at the Achziv Canyon, 32°42. 2631' 
North, 34°33. 5556' East at 1073 m. A subsample from Dor 
disturbance sediment was frozen for molecular analysis. All other 
samples were collected during the 201 1 Nautilm E/V field season. 
We used 70 mm diameter, 30 cm length cores to sample 9-10 cm 
profiles at the Dor disturbance location near shore of Israel. A core 
was taken from the yellow biofilm patch at 32°35. 3227' North, 
34°42. 9271' East at a depth of 567 m. The biofilm was coUected 
with a pipette and flash-frozen and its pore-water was analyzed for 
Fe 2+ . Additionally, a biofilm and a control core were collected at 
32°34. 7468' North, 34°44. 4832' East at a depth of 379 m. These 
cores were sliced to 1 cm sections, and the sediment was flash- 
frozen for further processing. 

Microscopy 

The flocculent matter from the 567 m core was transferred with 
a pipette to a slide and observed immediately under various 
magnifications. The structure of the biofilm was observed with an 



Imager M2 (Carl Zeiss, Germany) microscope. The images were 
taken with an AxioCam MRm camera (Carl Zeiss, Germany). 

pH and 0 2 measurements 

The pH and 0 2 measurements were performed immediately 
upon the retrieval of the 567 m core. The intact core with the in 
situ headspace water was subjected to a gentle air bubbling to 
maintain steady state oxygenation. The FeOB biofilms were easily 
disturbed during manipulation; hence the profiles were not taken 
specifically through the biofilm. The pH and 0 2 concentration 
within the sediment were measured with pH-100 and O 2 -100 
microsensors (Unisense, Denmark), respectively, according to 
manufacturer's instructions. We utilized a manually operated 
micromanipulator to obtain the profile at 1 mm resolution. The 
pH and 0 2 concentrations were recorded at each depth following 
the stabilization of the signal. The signal was amplified with a 
Microsensor Multimeter (Unisense, Denmark) and processed in 
the SensorTrace Basic 3.0 software (Unisense, Denmark). Only 
one profile was obtained to allow immediate further processing of 
the core. 

Fe 2+ and total Fe determination 

The pore fluids were extracted by centrifugation under N 2 
atmosphere. Filtered subsamples were assayed immediately by 
ferrozine reagent [30] and Fe 2+ in the pore water was determined 
spectrophotometrically [31]. Total Fe in the biofilm was also 
determined with the ferrozine assay [31] as follows: 100 ul biofilm 
subsample was acidified with 1 mol 1 _1 HC1 overnight, followed 
by reduction with hydroxylamine hydrochloride prepared in a 
solution of analytical grade HC1. Following the ferrozine reagent 
addition, ammonium acetate solution adjusted to pH 9.5 was 
added to buffer pH. Only technical replicates were performed, due 
to limited availability of samples. The sensitivity of this assay is 
0.1 nmol l -1 [32], hence our results are several orders of 
magnitude above the detection limit. 

DNA isolation, amplification and sequencing 

DNA was isolated from flash-frozen sediments or biofilm, using 
PowerSoil DNA Isolation Kit (MoBio, Carlsbad, CA, USA) 
following manufacturer's instructions. The template DNA was 
amplified using GoTaq Green Master Mix (Promega, Madison, 
WI, USA) with addition of 1 ju.1 bovine serum albumin per 50 |0.1 
reaction in T100 thermal cycler (Bio-Rad, Hercules, CA, USA). 
The 16S rRNA gene was amplified using 27F-CM and 1492-R 
primers [33], yielding — 1 450 bp product; the 23S rRNA gene was 
amplified with broad specificity forward primer L-0858-a-S-21 
and zetaproteobacteria specific L-C-Zeta-161 l-A-22 reverse 
primer [6], yielding ~750 bp product. The PCR conditions were 
as follows: denaturation at 94°C for 30 s, annealing at for 30 s 50 
°C and extension at 72°C for lmin (1.5 min for 16S) for 35 cycles. 
PCR products were purified using the Wizard SV Gel and PCR 
Clean-Up System (Promega, Madison, WI, USA) and subsequent- 
ly cloned in the pGEM-T Easy vector (Promega, Madison, WI, 
USA) using Escherichia coli JM109 (Promega, Madison, WI, USA) 
as a host. The sequencing was performed by HyLabs (Israel) using 
pGEM-T specific primers T7 and SP6. 

Bacterial tag-encoded FLX amplicon pyrosequencing 
(bTEFAP) and data analysis 

16S rRNA gene was amplified using 16S universal Eubacterial 
primers 104F: 5' GGC GVA CGG GTG AGT AA; and 530R: 5' 
CCG CNG CNG CTG GCA C [34]. A single-step 30 cycle PCR 
using HotStarTaq Plus Master Mix Kit (Qiagen, Valencia, CA, 
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Figure 1. The map of the study sites, a) The bathymetric chart of the Levantine basin displaying areas explored in this study. Yellow circles 
indicate deposit covered sediments. Yellow to black gradient indicates presence of both orange and black patches. Black circles indicate presence of 
only black patches. White circles mark areas where no deposits were observed, b) A detailed track of the Dor disturbance survey. The sediment 
sample locations are shown. 
doi:10.1371/journal.pone.0091456.g001 



USA) were used under following conditions: denaturation at 94 °C 
for 30 s; annealing at 53 °C for 40 s and elongation at 72 °C for 1 
min for 28 cycles. Following PCR, all amplicon products from 



different samples were mixed in equal concentrations and purified 
using Agencourt Ampure beads (Agencourt Bioscience Corpora- 
tion, MA, USA). The samples were sequenced utilizing 454 GS 
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FLX titanium (Roche, Penzberg, Germany) following manufac- 
turer's guidelines. The data derived from the sequencing was 
processed using a proprietary analysis pipeline [35-38] at MR 
DNA (Shallowater, TX, USA). Barcodes and primers were deleted 
from the sequences, short sequences < 200bp were removed, 
sequences with ambiguous base calls were removed, and sequences 
with homopolymer runs exceeding 6bp were removed. Sequences 
were denoised and chimeras were removed using custom software 
[39]. Operational taxonomic units were defined after removal of 
singleton sequences, clustering at 3% divergence (97% similarity). 
OTUs were then taxonomically classified using BLASTn against a 
curated Greengenes database [40] and compiled into each 
taxonomic level. A total of 1659 valid sequences were established, 
yielding 97 OTUs. 

Phylogenetic analysis 

Evolutionary history was deduced using the maximum likeli- 
hood method based on Kimura 2-parameter model [41] (+/ for 
454 16S rRNA genesequences (333 positions), +G for 23S rRNA 
gene sequences (798 positions)) and Data Specific model [42] (+/, 
1298 positions) for the cloned 16S rRNA gene sequences. The 
bootstrap consensus tree inferred from 1000 replicates is taken to 
represent the evolutionary history of the taxa analyzed [43] . Initial 
trees for the heuristic search were obtained automatically as 
follows: When the number of common sites was less than 100 or 
less than one-fourth of the total number of sites, the maximum 
parsimony method was employed; otherwise BIONJ method with 
MCL distance matrix was used. Codon positions included were 
1 st+2nd+3rd+Noncoding. Evolutionary analyses were conducted 
in MEGA5 [44]. All sequences were submitted to GeneBank 
under accession numbers KF199322-KF199336, KF651136- 
KF651143. 

Results 

The distribution of biofilm patches deduced from the 
ROV imagery 

Sediments at 300-800 m areas at the Achziv canyon, Dor 
disturbance and the Palmachim disturbance (Fig. 1), were covered 
by sporadically distributed patches of yellow-orange colored 
flocculent material (Fig. 2a). The size of each biofilm varied from 
1 to 5 cm, yet detached particles were scattered over an 
approximated 10 cm diameter from the biofilm center (Fig. 2b). 
This phenomenon was consistent throughout the ca. 5 km surveys 
at all explored locations of the latter depth range. The maximum 
density of the patches was 7 units m 2 , resulting in approximately 
3% coverage of the sediment surface by the biofilm. We noticed 
that the biofilm was easily disturbed and moved by the ROV's 
motion or by biological activity, such as fish movement (Fig. 2c, 
supplementary video). Moreover, we observed the loss of 
coloration at depths greater than 800 m at the Achziv canyon, 
Acre and Palmachim locations (Fig. 2d). Below 800 m the 
coloration of patches turned black, the biofilm material became 
less flocculent and no longer easily disturbed (see Fig. S2 for the 
micrograph of the latter). 

Microscopic observation 

At 8x magnification, we observed packed matter with rusty 
coloration (Fig. 2f). A complex matrix was revealed after a closer 
inspection at 63x magnification, including transparent organic 
filaments and darker matrix of higher density (Fig. 2g, h). Small 
(less than 1 u.m) particles were attached along the matrix. Some 
details from this matrix resembled the stalks with attached iron 




Figure 2. Images of the FeOB patchy mats observed in this 
study, a) Biofilm patches at the sediment-water interface of Dor 
disturbance, b) In situ close up image of the deposits, the distance 
between two red dots is 10 cm. c) Dipturus sp. re-suspending the 
deposits, d) Black-colored deposits in at 1000 m depth offshore Israel, e) 
Close up image of the FeOB patch. The arrow points at potential 
burrow opening at the core of the patch, f) Microscopic image of the 
deposits ex situ (scale bar is 0.5 mm), g, h) Microscopic images of the 
deposits (scale bar is 20 |im). 
doi:10.1371/journal.pone.0091456.g002 

oxide particles, as found in the Fe-oxidizing zetaproteobacteria 
[6,22,45]. 

Physical and chemical parameters of the underlying 
sediments 

The oxygen profile was determined only at the top 6 mm layer 
in the intact core retrieved from 567 m Dor disturbance sediment. 
The oxygen dropped from 187 umol L 1 to 83 uxnol L . The 
pH also decreased throughout the 6 mm profile, from 7.72 to 
7.40, stabilized from 4 to 6 mm (Fig. 3a). Fe 2+ concentration was 
0.2-0.4 u.mol L 1 at first 2 cm of the sediment, followed by a 
maximum of 6.7 umol L 1 at 2.5 cm below sediment-water 
interface (Fig. 3b). The total acid leachable iron concentration 
within the biofilm reached a value of 15.8 ±0.1 mmol L _1 . 

Zetaproteobacteria within the biofilm at the sediment- 
water interface 

In order to establish the phylogenetic relationship of the main 
biofilm bacteria, ca. 1450 bp 16S rRNA gene fragments derived 
from the Sanger sequencing were analyzed. This allowed for more 
accurate phylogenetic clustering of the studied bacteria than a 
phylogenetic analysis of shorter (ca. 400 bp) sequences derived 
from bTEFAP [46] . The analysis of the bacterial population from 
the biofilm revealed the dominance (9 out of 15 sequences) of 
zetaproteobacterial candidate, related (98% similarity) to the iron 
oxidizers found in back-arc hydrofhermal fields of the Southern 
Mariana Trough [5] (Fig. SI a). Moreover, it was identical to the 
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Figure 3. Geochemical parameters of the sampled sediments, a) pH and 0 2 profiles next to the sediment-water interface in the core collected 
at the depth of 567 m in Dor disturbance, b) Fe 2+ profile in the same core. (pH and 0 2 were measured at 1 mm resolution, Fe 2+ was measured at 1 cm 
resolution.) 

doi:1 0.1 371 /journal.pone.0091 456.g003 



zetaproteobacterial candidate sequence retrieved from a FeOB 
biofilm from the Dor disturbance in 2010. DNA extracted from 
surface biofilm was also amplified with zetaproteobacteria specific 
primers for 23S rRNA gene, resulting in determination and 
identification of this sequence in our sample (Fig. Sib), related to 
Mariprofundus sp. GSB2 (95% similarity) [6]. 

Biofilm bacterial population by bTEFAP 

Zetaproteobacterial PYROTAGSs dominated the biofilm 
(75.2% out of 1659 verified pyrotags) (Fig. 4). Two major 
zetaproteobacterial groups were derived based on clustering (Fig. 
5). One group, including 33% of the zetaproteobacterial pyrotags, 
clustered with Mariprofundus spp., while the second group, 
including 67% of the zetaproteobacterial pyrotags, clustered with 
sequences from hydrothermal fields of the Southern Mariana 
Trough [5] . Other common bacteria, representing more than 1 % 
of total bacterial sequences, included Ralstonia pickettii (Betaproteo- 
bacteria) - like pyrotags; three deltaproteobacterial pyrotags (one 
of them was most closely related to Magnetovibrio blakemorei); 
Lactococcus - like pyrotags; candidate division WS3 pyrotags and 
Cyanobacteria - like OTUs. Moreover, we have analyzed the 
relative abundance of bacterial pyrotags in sediment sections 1 cm 
below the biofilm and 1 cm sediment-water interface from the 
control core. No zetaproteobacterial pyrotags were detected in 
either of these sections, while alpha, delta and gammaproteobac- 
terial pyrotags were dominant, summing up to ~75% of total 
pyrotags in the respective samples (Fig. 4b, c). 

Discussion 

Our survey of Israel's continental margin benthos revealed that 
numerous yellow-orange patches existed within the 300-800 m 
depth range of the examined areas. These patches were defined by 
high total iron concentration and by the presence of iron-oxidizing 
zetaproteobacteria, previously described from hydrothermal vents, 
exposed oceanic crusts, hydrocarbon seeps and iron deposits. One 
study determined that zetaproteobacteria were present within 
brackish sediments [11]. This is the first report describing FeOB 
populations at deep continental margins lacking hydrothermal or 



hydrocarbon seepage. Although we were unable to determine 
conclusively the autotrophic nature of zetaproteobacterial biofilms, 
we suggest that these biofilms may provide added production to 
the benthic ecosystem based on the previously determined 
potential of zetaproteobacteria to fix CO z [1,2]. 

The phylogenetic analysis revealed two major clusters of 
zetaproteobacteria, one related to Mariprofundus genus, and a 
second more abundant cluster related to hydrothermal vent 
zetaproteobacteria from the Pacific Ocean. The metabolic 
preferences and factors dictating the relative abundance of the 
latter groups are still unclear. Other bacterial OTUs had relative 
abundance less than 2%, although the phylogenetic affiliation of 
some of them hints at the existence of iron-related microcosm. We 
have found pyrotags genetically identical to Ralstonia pickettii, 
bacterium previously reported from environments with high metal 
concentrations [47]. Moreover, pyrotags associated with magne- 
totactic genus Magnetovibrio [48] were found. Although the 
presence of cyanobacteria related pyrotags within the biofilm is 
surprising, their source can be the refractory part of marine snow. 
Moreover, several pyrotags were related to algal chloroplasts. We 
have previously determined similar sequences from the sediment- 
water interface at various locations in the deep eastern Mediter- 
ranean (data not shown). Hence, they may represent a large array 
of phytoplankton organisms, especially those protected by a hard 
cell wall due to an increased chance of preservation in the water 
column and higher sinking velocities [49] . 

Based on the hypothesis stating that the bioturbating/bioirrigat- 
ing organisms' burrows can induce the ferrous iron transport to 
the sediment surface [12], we can speculate that such burrows 
reach the sediment surface below the FeOB biofilm, providing 
reduced iron to zetaproteobacteria. Tubular hard structures, 
observed within the inactive black spot bulk sediment sample from 
1076 m at the Achziv Canyon (Fig. S2), are potential relict 
burrows with iron-manganese precipitate walls. This possibly 
explains the patchy appearance of the biofilm on the seafloor. The 
Fe 2+ is most likely produced in 2.5 cm below the sediment-water 
interface, where its concentration is ~ 30 times higher compared to 
0.5 cm below the biofilm, and its aided transport is a prerequisite 
for sufficient flux of Fe 2+ to sediment surface. 
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Figure 4. Relative abundance of major bacterial groups derived from bTEFAP analysis, a) Relative abundance of major microbial OTUs 
within the deposit collected from the 567 m Dor core, b) Relative abundance of major microbial classes within the 1 cm below the sediment-water 
interface section in the core taken from the deposit at the depth of 379 m in Dor disturbance, c) Relative abundance of major microbial classes within 
the 1 cm below the sediment-water interface section in the core taken from the non-deposit (control) sediment at the depth of 379 m in Dor 
disturbance. 

doi:1 0.1 371/journal.pone.0091 456.g004 
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Interestingly, at water depth of 800-1300 m biofilm patches 
were observed, although their black coloration and lack of 
flocculence may indicate loss of zetaproteobacterial activity and 
precipitation of manganese with iron oxides or pyrite formation. 
Several patches in transitional state (yellow to black) were also 
present (Fig. S3a, b). No biolilms were observed between 1300— 
1 700 m. The flux of organic matter and benthic infaunal biomass 
decrease with depth [50-52], therefore the increased organic load 



in shallow sediments can provide a larger pool of electrons for 
respiration compared to deeper sediments, producing larger 
quantities of reduced metabolites. This organic load can also 
support larger populations of bioturbating metazoans that may 
include sedentary Polychaeta families such as Spionidae, Aphar- 
etidae, Pectinariidae and Maldanidae, potentially responsible for 
the enhancement of Fe 2+ transport. Hence, a barrier disabling the 
sediments from providing reduced metabolites to the overlaying 
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EU491252.i l Uncultured bacterium clone P9X2b2C04 
JF973638.ll Uncultured bacterium clone MEBA6U1 
931 F3 17957. 1\ Mariprofundus sp. M34 

EF493243. 1 \ Mariprofundus ferrooxydans strain PV-1 

Israel Zela proteobaclerium pyrotag I (1 1) KF 199329 
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Hydrothermal group 



Zeta proteobacteria 
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Figure 5. Molecular phylogenetic analysis by maximum likelihood method of most common 16S rRNA pyrotags in the FeOB sample 
(grey color) and related sequences from the NCBI database. The number of pyrotags yielded for the respective OTU is shown in the brackets. 
Escherichia coli 16S sequence is used as an out-group. The numbers at nodes are bootstrap percent values based on 1000 resamplings. The scale bar 
corresponds to the number of substitutions per site. 
doi:1 0.1 371 /journal.pone.0091 456.g005 
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water can be produced by insufficient organic matter flux to the 
seafloor. The seasonal fluctuations dependent on the primary 
productivity in the photic zone may explain the inactive iron 
oxidation state of 800-1300 m sediments. Both 2010 and 2011 
legs of the Nautilus E/V exploration in Levantine basin took place 
in late summer prior to water column mixing and it should be 
noted that the depth of biofilm activity would be expected to 
increase following the spring phytoplankton bloom. 

Bioturbation/bioirrigation, sediment resuspension and chemical 
processes at the sediment-water interface were suggested previ- 
ously as important mechanisms dictating the fluxes of sediment 
iron into the water column [53]. In this study we established that 
the reduced Fe can reach the sediment surface and given the 
biological resuspension of sediment (Fig 2c, Supp. Video) and deep 
water mixing, can be transported up the water column. A 
prerequisite for the ability of iron to travel up the water column is 
the neutral buoyancy of the iron-containing particle. This can be 
achieved by detachment of small, colloidal iron particles. 
Moreover, the association of larger iron particles with positively 
buoyant organic matter can yield the same result [54,55]. Based 
on existing data [56] we estimate that winter sea surface 
temperatures (SST) below 15°C are sufficient to mix the water 
column to depths greater than 400 m. The average winter SST in 
the eastern Mediterranean fluctuates around 15.5 °C, while the 
lowest average temperature of 15°C was detected during the 
winter of 1992 [56]. During the last decade, there is an increasing 
trend of SST values in the eastern Mediterranean Sea [56] . The 
enhanced SST is expected to decrease iron supply from the 
sediment surface due to shallower mixing depth [57]. On the 
contrary, increased temperature of eastern Mediterranean deep 
water sources can result in the uplift of deep water [58] and 
increased nutrient flux from the sediment surface to the photic 
zone. Another positive effect of the global change on the nutrient 
supply may be the rare, extremely cold winter temperatures that 
were recently detected in this region [59]. Based on several 
models, their frequency may increase in the future [60]. In such 
case, the significance of the sediment-water interface iron as a 
bioavailable iron source for the water column may increase due to 
enhanced mixing during an extremely cold event. 

Unfortunately, water column iron profiles are still unavailable in 
the Levantine basin, and present Mediterranean sampling is 
limited to shallow (<10 m) or deep (>1500 m) stations [61,62]. 
Nonetheless, data from a Weddell Sea station, defined by 
proximity to the continental shelf and bathymetry resembling 
our study, suggests iron enrichment close to the sediment-water 
interface [63], verifying that sediment is a potential source for the 
water column iron. The bioavailability of the resuspended 
particulate iron originating from the sediment-water interface 
can be altered by biologically mediated processes such as microbial 
reductive dissolution [16,64-67] in the water column (Fig. 6). Iron 
minerals produced by the stalk forming FeOB are attached to 
carboxyl-rich polysaccharides [22], providing the organic ligands 
necessary for the reductive dissolution. 

In the report we show that iron-oxidizing bacteria precipitate 
the ferrous iron at the sediment-water interface along large areas 
of the continental margin, potentially altering its solubility. The 
fate of iron precipitated by the zetaproteobacteria is still unclear 
and its bioavailability may play a crucial role in the future of east 
Mediterranean ecology under the global change. 

Conclusions 

The "plain" sediments at continental margins constitute a 
heavily understudied area of deep sea research. In this study we 



Upper water column 




Mixing 



Reductive dissolution ? 




Figure 6. A conceptual model describing formation of the 
FeOB deposits and potential of the latter to contribute new 
iron to the water column. Orange to blue color gradient describes 
reductive process; blue to orange color gradient describes oxidative 
process. Briefly, Fe (III) oxides are important electron acceptors within 
the anoxic zone, forming Fe 2+ . Fe 2+ is diffusively transported to the 
upper oxygenated sediments, where it is being oxidized back to Fe 3+ . 
The bioturbating/bioirrigating organisms enhance the transport of Fe 2+ 
and allow it to reach the sediment surface, where it is oxidized by 
zetaproteobacterial FeOB. The flocculent FeOB deposits are easily 
disturbed and resuspended by larger fauna, being a source of 
particulate Fe to the water column. Particulate Fe can be reduced 
inside reducing microenvironments in the water column, becoming 
bioavailable. 

doi:10.1371/journal.pone.0091456.g006 

report the presence of patchy FeOB mats, abundant in the 
surveyed locations, unaffected by hydrothermal or hydrocarbon 
seepage, along the Levantine continental margin. Zetaproteobac- 
terial FeOB, represented by two main phylotypic groups derived 
from the 16S rRNA gene phylogenetic analysis, constitute the 
dominant fraction of bacterial population within the biofilms. We 
suggest that Fe 2+ concentration sufficient to maintain sediment- 
water interface FeOB community is most likely achieved by 
enhancement of pore water upward flux by bioturbating/ 
bioirrigating infauna. Given the potential of FeOB to fix carbon 
and alter iron chemistry, they may play an important role in the 
marine environment, affecting both benthic and water column 
ecosystems. 

Supporting Information 

Figure SI Molecular phylogenetic analysis by maxi- 
mum likelihood method of cloned bacterial genetic 
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markers, a) Evolutionary history of the 16S rRNA small subunit 
gene. Gamma-proteobacterial 16S sequence is used as an out- 
group, b) The evolutionary history of 23S rRNA large subunit 
gene. Mtrosomonas sp. 23S rRNA gene sequence is used as an out- 
group. The numbers are bootstrap percent values based on 1000 
resamplings. The scale bar corresponds to the number of 
substitutions per site. 
(TIF) 

Figure S2 Microscopic images of precipitated material 
within sediments from a black spot, a) Biogenic tubes, b) 
Solid particles. Scale bar is 0.5 mm. 
(TIF) 

Figure S3 Additional images of FeOB deposits, a) A 

deposit in transitional state, b). Deposits next to a metal wreck 
have both yellow-orange and black colorations, c) A spot in the 
moddle of FeOB deposit that can be an opening of metazoan 
burrow, d) Disturbance and resuspention of FeOB deposit by 
deep-sea Pleuronectiformes specimen. 
(TIF) 

Video SI Re-suspension of deposit-covered sediment by 
Dipturus specimen. 

(MOV) 

References 

1. Emerson D, RcntzJA, Lilburn TG, Davis RE, Aldrich H, et al. (2007) A novel 
lineage of proteobaeteria involved in formation of marine Fc-oxidizing mierobial 
mat communities. PLoS One 2: c667. Available: http://www.plosonc.org/ 
articlc/info%3Adoi%2F10.1371%2Fjournal.pone.0000667. Accessed 2012 
March 9. 

2. Singer E, Emerson D, Webb EA, Barco RA, Kucncn JG, ct al. (2011) 
Mariprofundus fenocxydans PV-1 the first genome of a marine Fe(II) oxidizing 
Zctaproteobactcrium. PLoS One 6: c25386. Available: http:/ /www.plosone. 
org/article/info% 3Adoi%2F 10.1371 %2Fjournal.pone.0025386. Accessed 20 1 2 
March 1. 

3. Dang H, Chen R, Wang L, Shao S, Dai L, ct al. (2011) Molecular 
characterization of putative biocorroding microbiota with a novel niche 
detection of Epsilon- and Zetaproteobaeteria in Pacific Ocean coastal seawaters. 
Environ Microbiol 13: 3059-3074. 

4. McAllister SM, Davis RE, McBeth JM, Tebo BM, Emerson D, ct al. (2011) 
Biodiversity and emerging biogeography of the neutrophilic iron-oxidizing 
Zetaproteobaeteria. Appl Environ Microbiol 77: 5445-5457. 

5. Kato S, Yanagawa K, Sunamura M, Takano Y, Ishibashi J, et al. (2009) 
Abundance of Zetaproteobaeteria within crustal lluids in back-are hydrothermal 
fields of the Southern Mariana Trough. Environ Microbiol 11: 3210-3222. 

6. McBeth JM, Little BJ, Ray RI, Farrar KM, Emerson D (2011) Neutrophilic 
iron-oxidizing "zetaproteobaeteria" and mild steel corrosion in nearshore 
marine environments. Appl Environ Microbiol 77: 1405-1412. 

7. Omorcgic EO, Mastalcrz V, dc Langc G, Straub KL, Kappler A, et al. (2008) 
Biogeochcmistry and community composition of iron- and sulfur-precipitating 
mierobial mats at the Chefren mud volcano (Nile Deep Sea Fan, Eastern 
Mediterranean). Appl Environ Microbiol 74: 3198-3215. 

8. Kato S, Nakamura K, Toki T, Ishibashi J-I, Tsunogai U, et al. (2012) Iron-based 
mierobial ecosystem on and below the sealloor: a ease study of hydrothermal 
fields of the southern mariana trough. Front Microbiol 3: 89. Available: http:// 
www.pubmedcentral.nili.gov/ articlerender.fcgi?artid — 3304087&tool — pmccntrcz& 
rendertype — abstract. Accessed 2012 Oct 19. 

9. Emerson D, Fleming EJ, McBeth JM (2010) Iron-oxidizing bacteria: an 
environmental and genomic perspective. Annu Rev Microbiol 64: 561—583. 

10. Emerson D, Moyer C (2010) Microbiology of seamounts: common patterns 
observed in community structure. Oceanography 23: 148-163. 

1 1 . McBeth JM, Fleming EJ, Emerson D (2013) The transition from freshwater to 
marine iron-oxidizing bacterial lineages along a salinity gradient on the 
Sheepscot River, Maine USA. Environ Microbiol Rep 5: 453—463. 

12. Scvcrmann S, McManus J, Bcrclson WM, Hammond DE (2010) The 
continental shelf benthic iron flux and its isotope composition. Geoehim 
Cosmochim Acta 74: 3984-4004. 

13. Sundby B (2006) Transient state diagencsis in continental margin muds. Mar 
Chem 102: 2-12. 

14. Chase Z (2005) Distribution and variability of iron input to Oregon coastal 
waters during the upwelling season. J Geophys Res 110: C10S12. Available: 
http://www.agu.org/pubs/crossref/2005/2004JC002590.shtml. Accessed 2012 
Nov 11. 

15. Elrod VA. (2004) The flux of iron from continental shelf sediments: A missing 
source for global budgets. Geophys Res Lett 31: L12307. Available: http:// 



Acknowledgments 

We thank the European Commission FP7 research infrastructure initiative 
program, "Assemble 227799," for partial support of this project. This 
research used samples and data provided by the E/V Nautilus Exploration 
Program - expeditions, NA009 and NA019. The authors would like to 
thank all individuals who helped during the expedition, including onboard 
technical and scientific personnel, and the captain and crew of the E/V 
Nautilus. We also thank the reviewers for constructive comments on an 
earlier version of this manuscript. Supplementary information is available 
at Plos ONE website. All authors of this manuscript certify that they qualify 
for authorship because of substantial contribution to the work submitted. 
The authors declare that this manuscript has not been published nor is 
under simultaneous consideration for publication elsewhere. The authors 
agree to transfer the copyright to the Plos ONE journal to be effective if 
and when the manuscript is accepted for publication and that the 
manuscript will not be published elsewhere in any other language without 
the consent of the Plos ONE journal. 

Author Contributions 

Conceived and designed the experiments: MRB GA RT BGTJA DC ZBA 
DT. Performed the experiments: MRB GA RT ES DT. Analyzed the data: 
MRB GA RT ES DT. Wrote the paper: MRB GA RT ES JA DC ZBA 
BGT DT. Saw and approved the final manuscript: MRB GA RT ES JA 
DC BGT ZBA DT. 



www.agu.org/pubs/crossrcf/2004/2004GL020216.shtml. Accessed 2012 Nov 
9. 

16. Lohan MC, Bruland KW (2008) Elevated Fe(II) and dissolved Fe in hypoxic 
shelf waters off Oregon and Washington: an enhanced source of iron to coastal 
upwelling regimes. Environ Sci Technol 42: 6462—6468. 

17. Lam PJ, Bishop JKB, Henning CC, Marcus MA, Waychunas GA, et al. (2006) 
Wintertime phytoplankton bloom in the subarctic Pacific supported by 
continental margin iron. Global Biogeochem Cycles 20: 1-12. 

18. Lam PJ, Bishop JKB (2008) The continental margin is a key source of iron to the 
HNLC North Pacific Ocean. Geophys Res Lett 35: L07608. Available: http:// 
www.agu.org/pubs/crossref/2008/2008GL033294.shtml. Accessed 2012 Nov 
11. 

19. Katz T, Yahcl G, Rcidenbach M, Tunniclifle V, Hcrut B, et al. (2012) 
Resuspension by fish facilitates the transport and redistribution of coastal 
sediments. Limnol Oceanogr 57: 945-958. 

20. Rich HW, Morel FMM (1990) Availability of well-defined iron colloids to the 
marine diatom Uidassiosim wdssflogh. Limnol Oceanogr 35: 652—662. 

21. Toner BM, Berquo TS, Michel FM, SorensenJ V, Templcton AS, ct al. (2012) 
Mineralogy of iron microbial mats from loihi seamount. Front Microbiol 3: 1 18. 
Available: http:/ /www .pubmedccntral.nih.gov/ artielerender .fcgi?artid = 3316996& 
tool = pmccntrcz&rendcrtype = abstract. Accessed 2013 Dec 17. 

22. Chan CS, Fakra SC, Emerson D, Fleming EJ, Edwards KJ (2011) Lithotrophic 
iron-oxidizing bacteria produce organic stalks to control mineral growth: 
implications for biosignature formation. ISME J 5: 717—727. 

23. Tanaka T, Zohary T, Krom MD, Law CS, Pitta P, ct al. (2007) Microbial 
community structure and function in the Levantine Basin of the eastern 
Mediterranean. Deep Sea Res Part I Oceanogr Res Pap 54: 1721—1743. 

24. Bcrman-Frank I, CuUenJT, Shakcd Y, ShcrrcU RM, Falkowski PG (2001) Iron 
availability, cellular iron quotas, and nitrogen fixation in Tridwdesmium. Limnol 
Oceanogr 46: 1249-1260. 

25. Brcitbarth E, Achterberg EP, Ardelan MV, Baker AR, Bucciarelli E, et al. (20 1 0) 
Iron biogeochcmistry across marine systems - progress from the past decade. 
Biogeosciences 7: 1075-1097. 

26. Yogcv T, Rahav E, Bar-Zecv E, Man-Aharonovich D, Stamblcr N, ct al. (201 1) 
Is dinitrogen fixation significant in the Levantine Basin, East Mediterranean 
Sea? Environ Mcrobiol 13: 854-871. 

27. Thcodosi C, Markaki Z, Mihalopoulos N (2010) Iron speciation, solubility and 
temporal variability in wet and dry deposition in the Eastern Mediterranean. 
Mar Chem 120: 100-107. 

28. Coleman D, Austin JrJA, Ben-Avraham Z (2011) Exploring the Continental 
Margin of Israel: " Telepresence " at Work. Eos Trans AGU 92. Available: 
http://www.agu.org/pubs/crossref/20 1 1 /20 1 1EO 1 00002. shtml. Accessed 
20 12 June 24. 

29. Coleman DF, Austin JrJA, Ben-Avraham Z, Makovsky Y, Tchcrnov D (2012) 
Seafioor pockmarks, deepwater corals, and cold seeps along the continental 
margin of Israel. Oceanography 25 (supplement 1): 40-41. 

30. Stookey LL (1970) Ferrozinc-a new spectrophotometrie reagent for iron. Anal 
Chem 42: 779-781. 



PLOS ONE | www.plosone.org 



9 



March 2014 | Volume 9 | Issue 3 | e91456 



Zetaproteobacterial Mats, Levant 



31. Viollicr E, Inglctt P, Hunter K, Roychoudhury A, Van Gappcllcn P (2000) The 
ferrozinc method revisited: Fe(II)/Fe(III) determination in natural waters. Appl 
Geochemistry 15: 785-790. 

32. Achterberg EP, Holland TW, Bowie AR, Mantoura RFC, Worsfold PJ (2001) 
Determination of iron in seawater. Anal Chim Acta 442: 1-14. 

33. Frank JA, Reich CI, Sharma S, Weisbaum JS, Wilson BA, et al. (2008) Critical 
evaluation of two primers commonly used for amplification of bacterial 16S 
rRNA genes. Appl Environ Microbiol 74: 2461-2470. 

34. Wang Y, Qian P (2009) Conservative fragments in bacterial 16S rRNA genes 
and primer design for 16S ribosomal DNA amplicons in metagenomic studies. 
PLoS One 4: c7401. Available: http://ww^v.plosonc.org/artiflc/info%3Adoi% 
2F10.1371%2Fjournal.pone.0007401. Accessed 19 August 2013. 

35. Capone KA, Dowd SE, Stamatas GN, Nikolovski J (2011) Diversity of the 
human skin microbiomc early in life. J Invest Dermatol 131: 2026—2032. 

36. Dowd SE, Sun Y, Wolcott RD, Domingo A, GarroUJA (2008) Bacterial tag- 
encoded FLX amplicon pyrosequencing (bTEFAP) for microbiomc studies: 
bacterial diversity in the ileum of newly weaned Salmonella-infected pigs. 
Foodborne Pathog Dis 5: 459-472. 

37. Eren AM, Zozaya M, Taylor CM, Dowd SE, Martin DH, et al. (201 1) Exploring 
the diversity of Gardnerella vaginalis in the genitourinary tract microbiota of 
monogamous couples through subtle nucleotide variation. PLoS One 6: e26732- 
Available: http://www.plosone.org/article/inib%3Adoi%2F10. 1371 %2Fjournal. 
ponc.0026732. Accessed 2012 Nov 12. 

38. Swanson KS, Dowd SE, Suchodolski JS, Middelbos IS, Vester BM, ct al. (201 1) 
Phylogcnctic and gene-centric mctagenomics of the canine intestinal micro- 
biome reveals similarities with humans and mice. ISME J 5: 639-649. 

39. Gontcharova V, Youn E, Sun Y, Wolcott RD, Dowd SE (20 1 0) A comparison of 
bacterial composition in diabetic ulcers and contralateral intact skin. Open 
MicrobiolJ 4: 8-19. 

40. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, et al. (2006) 
Greengenes, a chimera-checked 16S rRNA gene database and workbench 
compatible with ARB. Appl Environ Microbiol 72: 5069-5072. 

41. Kimura M (1980) A simple method for estimating evolutionary rate of base 
substitutions through comparative studies of nucleotide sequences. J Mol Evol 
16: 111-120. 

42. Nei M, Kumar S (2000) Molecular Evolution and Phylogenctics. New York: 
Oxford University Press. 

43. Felscnstein J (1985) Confidence limits on phylogcnics: An approach using the 
bootstrap. Evolution (N Y) 39: 783-791. 

44. Tamura K, Peterson D, Peterson N, Stechcr G, Nei M, ct al. (201 1) MEGA5: 
molecular evolutionary genetics analysis using maximum likelihood, evolution- 
ary distance, and maximum parsimony methods. Mol Biol Evol 28: 2731—2739. 

45. Hodges TW, Olson JB (2009) Molecular comparison of bacterial communities 
within iron-containing flocculent mats associated with submarine volcanoes 
along the Kermadcc Arc. Appl Environ Microbiol 75: 1650—1657. 

46. Wommack KE, Bhavsar J, Ravel J (2008) Mctagenomics: read length matters. 
Appl Environ Microbiol 74: 1453-1463. 

47. Konstantinidis KT, Isaacs N, Fctt J, Simpson S, Long DT, et al. (2003) 
Microbial diversity and resistance to copper in metal-contaminated lake 
sediment. Microb Ecol 45: 191-202. 

48. Bazylinski DA, Williams TJ, Lefevrc CT, Trubitsyn D, Fang J, et al. (2013) 
Magnetovikrio klakemorei gen. nov., sp. nov., a magnetotactic bacterium 
(Alphaproteobacteria: Rhodospirillaccae) isolated from a salt marsh. Int J Syst 
Evol Microbiol 63: 1824-1833. 



49. Iversen MH, Ploug H (2010) Ballast minerals and the sinking carbon flux in the 
ocean: carbon-specific respiration rates and sinking velocity of marine snow 
aggregates. Biogeoscicnccs 7: 2613-2624. 

50. Rijk S De, Jorissen FJ, Rohling EJ, Troelstra SR (2000) Organic flux control on 
bathymctric zonation of Mediterranean benthic foraminifcra. Mar Micropa- 
lcontol 40: 151-166. 

51. Danovaro R, Croce N Delia, Eleftheriou A, Fabiano M, Papadopoulou N, et al. 
(1995) Meiofauna of the deep Eastern Mediterranean Sea: distribution and 
abundance in relation to bacterial biomass, organic matter composition and 
other environmental factors. Prog Oceanogr 36: 329-341. 

52. Danovaro R, Dinet A, Duineveld G, Tselcpides A (1999) Benthic response to 
particulate fluxes in different trophic environments: a comparison between the 
Gulf of Lions-Catalan Sea (western-Mediterranean) and the Cretan Sea 
(eastern-Mediterranean). Prog Oceanogr 44: 287—312. 

53. Pakhomova SV, Hall POJ, Kononets MY, Rozanov AG, Tengberg A, et al. 
(2007) Fluxes of iron and manganese across the sediment-water interface under 
various redox conditions. Mar Chcm 107: 319—331. 

54. Azetsu-Scott K, Passow U (2004) Ascending marine particles: Significance of 
transparent exopolymer particles (TEP) in the upper ocean. Limnol Oceanogr 
49: 741-748. 

55. Ortega-Osorio A, Scott SD (2001) Morphological and chemical characterization 
of neutrally buoyant plume-derived particles at the Eastern Manus basin 
hydrothermal field, Papua New Guinea. Can Mineral 39: 17—31. 

56. Hcrut B, Almogi-Labin A, Jannik N, Gcrtman I (2000) The seasonal dynamics of 
nutrient and chlorophyll a concentrations on the SE Mediterranean shelf-slope. 
Occanol Acta 23: 771-782. 

57. Coma R, Ribes M, Serrano E, Jimenez E, Salat J, et al. (2009) Global warming- 
enhanced stratification and mass mortality events in the Mediterranean. Proc 
Natl Acad Sci 106: 6176-6181. 

58. Danovaro R, Dell' Anno A, Fabiano M, Pusceddu A, Tselepides A (2001) Deep- 
sea ecosystem response to climate changes: the eastern Mediterranean case 
study. Trends Ecol Evol 16: 505-510. 

59. Tolika K, Anagnostopoulou C, Mahcras P, Velikou K (2013) Extreme 
temperature contrast of the year 2012 in Greece: An exceptionally cold winter 
and a record breaking summer. EGU General Assembly 2013. p. 2922. 

60. Sanchez E, Gallardo C, Gacrtner MA, Arribas A, Castro M (2004) Future 
climate extreme events in the Mediterranean simulated by a regional climate 
model: a first approach. Glob Planet Change 44: 163-180. 

61. Guieu G (2002) Impact of high Saharan dust inputs on dissolved iron 
concentrations in the Mediterranean Sea. Gcophys Res Lett 29: 191 1. Available: 
http://doi.wiley.com/10.1029/2001GL014454. Accessed 2013 June 3. 

62. Sarthou G, Jeandel C (2001) Seasonal variations of iron concentrations in the 
Ligurian Sea and iron budget in the Western Mediterranean Sea. Mar Chem 74: 
115-129. 

63. Westerlund S, Ohman P (1991) Iron in the water column of the Weddell Sea. 
Mar Chem 35: 199-217. 

64. Shaked Y, Lis H (2012) Disassembling iron availability to phytoplankton. Front 
Microbiol 3: 123. 

65. Boyd PW, Ellwood MJ (2010) The biogeochemical cycle of iron in the ocean. 
Nat Geosci 3: 675-682. 

66. Bruland K, Orians K, Cowcn J (1994) Reactive trace metals in the stratified 
central North Pacific. Geochim Cosmochim Acta 58: 3171—3182. 

67. Wu J, Boyle E (2002) Iron in the Sargasso Sea: Implications for the processes 
controlling dissolved Fe distribution in the ocean. Global Biogeochcm Cycles 16: 
33-1-33-38. 



PLOS ONE | www.plosone.org 



10 



March 2014 | Volume 9 | Issue 3 | e91456 



